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Abstract. Deeper understanding of the properties of dark energy via SNIa surveys, 
and to a large extent other methods as well, will require unprecedented photometric pre- 
cision. Laboratory and solar photometry and radiometry regularly achieve precisions 
on the order of parts in ten thousand, but photometric calibration for non- solar astron- 
omy presently remains stuck at the percent or greater level. We discuss our project to 
erase this discrepancy, and our steps toward achieving laboratory-level photometric pre- 
cision for surveys late this decade. In particular, we show near-field observations of the 
balloon-borne light source we are presently testing, in addition to previous work with 
a calibrated laser source presently in low-Earth orbit. Our technique is additionally ap- 
plicable to microwave astronomy. Observation of gravitational waves in the polarized 
CMB will similarly require unprecedented polarimetric and radiometric precision, and 
we briefly discuss our plans for a calibrated microwave source above the atmosphere as 
well. 



1. INTRODUCTION 

The presence of artificial flux standards above the Earth's atmosphere may provide sig- 
nificant reduction of photometric uncertainties for measurements that depend on such 
calibration. The combined effect of atmospheric and instrumental extinction in the vis- 
ible and near-infrared is presently the source of the largest uncertainty on the properties 
and amount of dark energy (see Fig. Q]). Man-made visible and NIR light sources can 
be measured to a precision of up to 100 times better than standard stellar sources, as 
presently shown by both laboratory and solar irradiance measurements, allowing for 
vast reduction of astronomical and cosmological uncertainties due to photometry. 

Separately, in the microwave spectrum, polarization patterns in the cosmic mi- 
crowave background (CMB) encode a vast amount of other information (beyond dark 
energy) about the early universe, ranging from the amplitude of primordial gravita- 
tional waves and the energy scale of cosmic inflation, to the gravitational lensing po- 
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Figure 1. Constraints on the dark energy equation of state parameters wo (= pip 
of the universe at the present day), and w a = dw/da (also at present, where a is 
the scale factor of the universe), for all current and upcoming SNIa dark energy 
projects (plus the Planck CMB mission) with the new calibration provided by our 
calibration program, vs. with all current means of calibration. The left plot shows 
the uncertainties if one artificially constraints the universe to be flat, and the right 
plot shows the uncertainties with flatness constraint removed. Dashed lines show 
68% and 95% confidence int ervals for each case. The calibration will improve the 
dark energy "figure of merit" (lAlbrecht et al.l2006|) by a factor of 2.4, with additional 
improvement beyond 2014. 



tential integrated over cosmic history. The interpretation of these signals is entirely 
contingent on an accurate calibration of polarized instrumental sensitivity. As yet, no 
well-calibrated, polarized, celestial microwave sources exist, and in the coming gener- 
ation of microwave telescopes, this deficit may become a limiting factor in our ability 
to measure and understand the polarized microwave sky. Current calibration solutions 
rely on relatively nearby ground-based sources, requiring refocusing and special low- 
sensitivity detectors to handle the near-field source and additional atmospheric loading. 
The cleanest and simplest calibration solution would solve both of these issues, by loft- 
ing a source into the far field, at distances of order tens of kilometres on large microwave 
telescopes. 

Prior to the second half of the 20th century, the only sources of light above the 
Earth's atmosphere were natural in origin: stars, and reflected light from planets, moons, 
comets, etc. Natural sources have of course served extremely well in astronomy: 
through understanding the physical processes governing stellar evolution, we are now 
able to fairly precisely understand the spectra of stars used as calibration sources [e.g. 
Bohlin (2000)]. Nevertheless, in all stars the vast bulk of material, and the thermonu- 
clear processes that themselves provide the light, lie beyond our sight below the surface 
of the star. Superb models of stellar structure are available, but uncertainties of many 
types always remain. 

Since the launching of the first high-altitude balloons and satellites, a separate class 
of potential light sources in space and in near-space has become available. Observable 
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light from most satellites is primarily due to direct solar reflection, or reflection from 
Earth's albedo. While providing a convenient method of observing satellites, this light 
is typically unsuitable for use as a calibrated light source due to large uncertainties in 
the reflectivity (and, to a lesser extent, the precise orientation and reflective area) of 
satellites' surfaces^] 

However, balloon- and satellite-based calibration sources for ground-based tele- 
scopes are by no means technically prohibitive. As an example, a standard household 
25-watt tungsten filament lightbulb (which typically have a temperature of the order of 
3000 K and usually produce approximately 1 watt of visible light between 390 and 780 
nm) which radiates light equally in all directions from a 700 km low Earth orbit has an 
equivalent brightness to a 12.5-magnitude star (in the AB system, although for this ap- 
proximate value the system makes little difference). In general, the apparent magnitude 
of an orbiting lamp at a typical incandescent temperature which radiates isotropically 
is approximately given by 



where P is the power of the lamp in watts, and h is the height of the orbit in kilometers. 
The systematic uncertainty on the radiance of an optimally-designed above-atmosphere 
lamp, where cost is no object, would be dominated by the precision of radiometric 
monitoring technology, and be in the range of approximately 200 parts per million if 
the best presently available radiometric technology is used. 

An alternative to an isotropic or near-isotropic lamp would be a laser source, with 
beam pointed at the observer (with a small moveable mirror, for example). Divergences 
of laser beams are typically on the order of a milliradian (which can be reduced to 
microradians with a beam expander) so much less output power than a lamp would 
be required for a laser beam to mimic the brightness of a typical star. The apparent 
magnitude of an orbiting laser with Gaussian beam divergence pointed directly at a 
ground-based telescope, is given by 



where P is the laser power in milliwatts, h is the height of the orbit in kilometers, and 
d is the RMS divergence of the laser beam in milliradians, under the assumption that 
the aperture of the telescope is small compared with the RMS width of the beam at the 
ground, hd. The RMS divergence would be the combination of the divergence at the 
source, and the divergence due to the atmosphere. In clear conditions, total atm ospheric 
diverg ence in a vertical path is at the level of approximately 5 microradians dTatarskil 
2006), and this of course only acts on the last fraction of the laser path that is within 
the atmosphere, so as long as the source divergence is significantly larger than this, at- 
mospheric divergence would be negligible. The uncertainty on the apparent magnitude 



Reflected solar light has, however, been successfully used as an absolute infrared calibration source by 
the Midcourse Space Experiment (MSX), using 2 cm diameter black-coated spheres ejected from th e MSX 
satellite, whose infrared emission was monitored by the instruments aboard MSX (Price et al. 2004]). This 
technique proved highly effective for the MSX infrared calibration; however, the technique is not easily 
applicable to measuring extinction of visible light in the atmosphere. 
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of an orbiting laser stemming from uncertainties in the radiometrically-monitored laser 
power would likely be limited by the precision of current radiometer technology. Mod- 
ern electrical substitution radiometers can achieve a precision of approximately 100 
parts per m illion when apertu re uncertainties can be neglected, as in the case of laser 
radiometry (iKopp et alJl2005h . Uncertainties on the magnitude due to uncertainty in the 



pointing and beam profile would potentially be limited by the size of the array of out- 
board telescopes for monitoring the laser spot, and by calibration differences between 
the individual telecopes in the array and with the main central telescope. The latter 
could clearly be minimized by a ground system for ensuring the relative calibration of 
the outboard telescopes and main telescope are all consistent. 

The uncertainties considered above assume that the exposure time is long com- 
pared with the coherence time of the atmosphere. With short exposures — or in the 
case of a laser that either quickly sweeps past, or is pulsed — atmospheric scintillation 
can play a major role in uncertainty in apparent magnitude of an above-atmosphere 
source. A typical timescale for a CW laser with 1 milliradian divergence in low Earth 
orbit to sweep past is tens of milliseconds, which is of the same order as characteristic 
timescales of atmospheric scintillation, and the typical timescale of single laser pulses 
is nanoseconds, much shorter than scintillation timescales, thus one cannot assume that 
such effects can be time-averaged over. In idealized conditions, for small apertures 
D <~ 5 cm and sub-millisecond integration times, the relative standard deviation in 
intensity 07 = AI/(I), where AI is the root-mean- square value of /, is given by the 
square root of 

oj = \9A2r 116 C 2 n (h)h 5/6 dh, (3) 
Jo 

where A is optical wavelength (in meters), Cl(h ) is known a s the refractive-index struc- 
ture coefficient, and h is altitude (in meters) (lTatarskill2006h . Large apertures D >~ 50 
cm have a relative standard deviation in intensity given by the square root of 

r*oo 

= 29.48ZT 7/3 \ C 2 n (h)h 2 dh (4) 



dTatarskil l2006h . The values and functional form of C 2 (h) are entirely dependent on 



the particular atmospheric conditions at the time of observation, however a relatively 
typical profile is given by the Hufnagel- Valley form: 

C 2 n (h) = 5.94 x lO^O^) 2 ^ V^ 1000 + 

2.7xlO- 16 e- h/15m +Ae- h/im , (5) 



where A and v are free parameters (Hufn agellll974l ). Commonly-used values for the 



A and v parameters, which represent the strength of turbulence near ground level and 
the high -altitude wind speed respectively, are A = 1.7 x 10" 14 m" 2/3 and v = 21 m/s 
(Rog gemann & W elsh 1996). Using these particular values, for a small aperture, the 



relative standard deviation 07 would be expected to be 0.466 for 532 nm light, which 
is n ot far off experi mental scintillation values for a clear night at a typical location 
[e.g. IJakemanl (Il978h l. For a single small camera, this is an extremely large uncertainty. 



Other than by increasing integration time (which is not possible with a pulsed laser) or 
by significantly increasing the camera aperture, the only way to reduce this uncertainty 
is to increase the number of cameras. With N cameras performing an observation, 
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which are spaced further apart than the coherence length of atmospheric turbulence 
(typically 5 to 50 cm), the uncertainty from scintillation can be reduced by a factor ^^N 
(for large N). 

The considerations above are necessarily both speculative and approximate. How- 
ever, at present there is an actual laser in low Earth orbit, visible with both equip- 
ment and with the naked eye, and analysis of ground-based observational data of the 
laser spot can be used for comparisons with the above, as well as for development of 
and predictions for potential future satellite-based photometric calibration sources of 
ground telescopes. 



2. CALIPSO Data and Analysis 



For c alibration of telescope optics and detector characteristics, authors (iStubbs et al 



2006) have both conceived of and used a wavelength-tunable laser within present and 



upcoming telescope domes as a color calibration standard. Further progress and details 
can be found within the proceedings from Stubbs & Tonry fr om this conference . Al- 
though a wavelength-tunable laser calibration source in orbit (I Albert et all 12006) does 
not exist yet, at present there is a 532 nm laser in low-Earth orbit pointed toward the 
Earth's surface, with precise radiometric measurement of the energy o f each of the 
20.25 Hz laser pulses, on the CALIPSO satellite, launched in April 2006 (iWinker et al 



2009b . We have collected data from a portable network of seven cameras and two cali- 



brated photodiodes, taken during CALIPSO flyovers on clear days in various locations 
in western North America. The cameras and photodiodes respectively capture images 
and pulses from the eye- visible green laser spot at the zenith during the moment of a 
flyover. Using precise pulse-by-pulse radiometry data from the CALIPSO satellite, we 
compare the pulse energy received on the ground with the pulse energy recorded by 
CALIPSO. The ratio determines the atmospheric extinction. 

The CALIPSO (Cloud Aerosol Lidar and Infrared Pathfinder Satellite Observa- 
tions) _satellite_^s_Jaunched on April 28, 2006 as a joint NASA and CNES mis- 
sion (iWinker et al.ll2009h . CALIPSO is part of a train of seven satellites (five of which 
are orbiting at the date of these proceedings), know n as the "A-Train," in sun - synchronous 
orbit at a mean altitude of approximately 690 km (iSavtchenko et al.ll2008h . CALIPSO 
completes an orbit every 98.4 minutes (approximately 14.6 orbits per day), and repeats 
its track every 16 days. CALIPSO contains a LIDAR (Light Detection and Ranging) 
system, known as CALIOP (Cloud Aerosol Lidar with Orthogonal Polarization), with 
a primary mission of obtain ing high resolutio n vertical profiles of clouds and aerosols 
in the Earth's atmosphere dHunt et al.ll2009h . The CALIOP laser produces simulta- 



neous, co-aligned 20 ns pulses of 532 nm and 1064 nm light, pointed a small angle 
(0.3°) away from the geodetic nadir in the forward along-track direction, at a repetition 
rate of 20.16 Hz. The light enters a beam expander, following which the divergence 
of each laser beam wavelength is approximately 100 yurad, producing a Gaussian spot 
of approximately 70 m RMS diameter on the ground. The pulse energy is monitored 
onboard the satellite, and averages approximately 110 mJ, at each one of the two wave- 
lengths, per pulse. The effective apparent magnitude of the 532 nm laser spot at the 
precise center of the beam is thus approximately -19.2, however this high brightness, 
of course, falls off rapidly as one moves away from the center of the beam. 

During 2007, the CALIPSO beam was observed at several locations in western 
North America using a portable ground station consisting of seven digital cameras, 
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Figure 3. Camera time-integrated irradiance data, and the resulting fitted time- 
integrated irradiance maps, for the (upper left) Apr. 17, (upper right) May 1, (lower 
left) May 30, and (lower right) Jun. 1 observations. The numbers at each camera 
refer to the measured value of time-integrated irradiance by that camera, with its 
associated uncertainty (68% CL), and the expectation value from the fitted function 
at the location of that camera. The contours on each plot are spaced at 1 /J/m 2 
intervals. The upper-right inset on each plot extends the x and y axes in order to see 
the three 2-D Gaussians that comprise the fitted function (as described in the text), 
and the lower-left inset shows a different 3-D view (from the side, rather than from 
above) of the fitted function and the data points. 
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and two calibrated photodiodes co-located with the central camera. Two example ob- 
servations can be seen in Fig. [2 Observing locations were selected by means of the 
CALIPSO ground track, monitored by NASA. 

The cameras were calibrated to an absolute radiometric standard through the use 
of the NIST-calibrated Hamamatsu photodiode (see above) and a low-noise 532 nm 
laser. Measurements of the response linearity, anisotropy, and ambient temperature 
dependence of the cameras were also performed. Linearity was found to be maintained 
to +2.4%. The anisotropy of the response of the CCDs was measured by obtaining 
images as above with the laser spot at 900 different places on each of the cameras' 
CCDs. The anisotropy, taken to be the standard devation of those 900 measurements, 
was found to be 0.5%. Temperature calibration was performed using one of the cameras 
inside a sealed refrigerator with a small hole in it to allow entrance of laser light (from 
inside a light-tight tube). Over a temperature range from 20° C down to 0°C, the image 
intensity was found to increase by the surprisingly large value of 52%. 

Knowing the precise relative positions of the cameras is necessary for extrapolat- 
ing the camera energy measurements into a measurement of the total laser pulse energy 
at the ground. During each observation, the location of each camera was measured us- 
ing GPS, as well as a surveyor's tape measure to determine the relative positions via 
triangulation. 

CALIPSO measures the energy of each of its individual laser pulses. The pulse en- 
ergy monitoring on CALIPSO consists of NIST-calibrated photodiodes mounted on an 
integrating sphere, with pulse-by-pulse energy measurement with design absolute pre- 



cision of +2% over full orbit, and relative precision of better than +0.4% (iWinker et al 



2009). The pulse energies are recorded in the CALIPSO datasets available from the At- 
mospheric Science Data Center (ASDC, located at NASA LaRCfl which are compared 
with the pulse energies measured at the ground to obtain measurements of atmospheric 
extinction. 

Plots of the resulting measurements can be found in Fig. [3l and the measured 
values of atmospheric extinction, limited in uncertainty by atmospheric extinctio n and 
the relatively poor knowledge of the shape of the laser beam, can be found in (lAlbert 
2012b . 



3. ALTAIR Balloon Program 

Beam-based pulsed sources are certainly not the only options for photometric calibra- 
tion standards above the atmosphere. More isotropic, and continuous, sources have the 
advantage that measurements are far less sensitive to the precise relative angle of the 
source and the observer, and that one can time-average over atmospheric scintillation, 
rather than having to aperture-average. Thus, a more isotropic calibrated source above 
the atmosphere is a more attractive option. No such sources presently exist. In order 
to test and fly such equipment, we have begun a program, denoted ALTAIR (Airborne 
Laser for Telescopic Atmospheric Interference Reduction), of flights and observations 
of a small high-altitude balloon calibrated source platform, shown in Fig.|U These bal- 
loon flights will take place over major surveys, such as Pan-STARRS and LSST, in 2014 
and beyond. Initial flights in New Hampshire have begun, and flights over Mt. Hopkins, 



2 http://eosweb.larc.nasa.gov 
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Figure 4. (Top row) 3-D solid diagram (left), photograph (center), and launch 
(right) of the balloon pay loads. (Second row) Setup for near- field imaging of the 
optical source: telescope (left) and view from above (right). (Third row) Near- field 
images of the optical source. (Bottom row) Commonly-available laser lines superim- 
posed on a plot of average atmospheric transmission (left) and measured light output 
distribution of the source (right). 
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Arizona will begin later this year. A balloon platform will not only provide critical data 
for precisely calibrating SNIa surveys, and also testing future satellite equipment, but 
additionally has the permanent advantage over satellite sources that the payload may 
be recovered following the flight, and tested in the laboratory, rather than only having 
such laboratory tests preceeding the launch. However, balloons of course cannot attain 
the altitude, nor the global reach, of a satellite. 

Our initial payload light source contains laser diode module sources at 440, 532, 
639, and 690 nm, directed into a 2" diameter integrated sphere via fiber optics, and 
monitored with a NIST-calibrated Hamamatsu S2281 photodiode. In close parallel with 
our program of test flights, we have also begun a campaign of near-field observations 
of the source, via placement of the source on a motorized mount on a rooftop, and 
observation with a telescope on a nearby rooftop approximately 75 m away, as shown 
in Fig. |U The addition of this man-made calibrated light source in near-space to the 
arsenal of techniques for photometric calibration will provide a powerful new tool for 
increasing precision in astrophysics. 

4. CONCLUSION 

We have performed initial tests and measurements of space- and near- space-based pre- 
cision photometric calibration sources. Our initial tests are highly promising, and we 
believe this is a viable means to reduce the dominant uncertainty in measurements 
of dark energy this decade. Our technique can additionally be applied to microwave 
astrophysics and to other regions of the spectrum, to impact cosmological and other 
measurements in those areas as well. Improved precision in photometric calibration 
will be nearly as critical for astronomy as increased aperture and etendue telescopes 
in upcoming decades. The future of precision photometry is extremely promising, and 
laboratory-based standards in near- space and in space allow one to for see many-fold 
improvement in photometric calibration as a near-term prospect. 
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